% stimulate a spike
% provided std, ts values
% hamming window forces the edges to be undisturbed
% latest date 22.4.3

function[ts] = spike(stddv,ts),

A = rand(1)*stddv; % standard deviation - idea of amplitude
W = hamming(length(ts)); % force ends to be near to values
ts = ts+((ts-mean(ts))*A).*W;


